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SYNCHROTRON-BETATRON  PARAMETRIC  INSTABILITY 
IN  FREE-ELECTRON  LASERS 


I.  Introduction 

All  electron  beams  have  finite  emittance;  that  is,  there  is  a  spread 
in  the  particle  velocities  transverse  to  the  primary  direction  of  motion. 

On  traversing  a  viggler,  the  inevitable  presence  of  transverse  gradients  in 
the  magnetic  field  lead  to  betatron  oscillations  of  the  particles.  On  the 
other  hand,  the  beating  of  the  viggler  and  optical  fields  forms  the 
ponderomotive  potential  wherein  the  particles  perform  synchrotron 
oscillations.  Coupling  of  these  two  basic  oscillatory  degrees  of  freedom 
in  a  free-electron  laser  (FEL)  may  be  deleterious  to  its  efficient 
operation  due  either  to  detrapping  out  of  the  ponderomotive  potential  well 
or  growth  of  betatron  oscillation  amplitude  and  emittance. 

In  Ref.  1  it  is  shown  that  the  curvature  of  the  optical  wavefronts 
couples  the  synchrotron  and  betatron  oscillations  of  the  electrons  and  that 
under  the  proper  resonance  condition  the  amplitude  of  the  synchrotron 
motion  increases  unstably.  In  subsequent  work,  the  authors  of  Ref.  2 
shoved  that  a  similar  betatron-induced  forcing  of  the  synchrotron  motion 
results  in  the  case  of  planar  optical  wavefronts  provided  the  viggler  field 
is  tapered. 

In  what  follows  it  is  shown  that  the  particle  motion  is  additionally 
susceptible  to  an  instability  which  is  dependent  neither  on  the  curvature 
of  the  wavefronts  nor  on  a  tapered  viggler  field.  The  instability  arises 
from  the  energy-dependence  of  the  betatron  wavenumber.  As  a  result  when  the 
electrons  undergo  synchrotron  oscillations  the  betatron  wavenumber  is 
modulated  at  the  synchrotron  period,  resulting  in  a  parametric  instability. 
This  leads  to  an  exponential  growth  of  the  betatron  oscillation  amplitude 
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and  velocity.  The  region  of  parameter  space  wherein  this  instability  is 
operative  is  identified  and  an  analytical  expression  for  the  growth  rate  is 
obtained. 

II.  Wiggle-average  equations  of  motion 

It  is  assumed  that  the  wiggler  is  linearly  polarized,  with  the  vector 
potential  given  by 

Aw  =  Awcosh(kwy)  sin(kwz)  ex, 

where  is  constant,  2n/kv  is  the  wiggler  period  and  ex  is  the  unit  vector 
along  the  x-axis.  The  vector  potential  of  the  linearly-polarized  optical 
field  is  taken  to  correspond  to  the  fundamental  vacuum  Gaussian  mode: 

Ag  =  Agcos(“  z  -  cat  +  <f>)  ex, 

where  «  is  the  angular  frequency,  c  is  the  speed  of  light  in  vacuo, 

♦(y.z)  =  «y2/(2cR)  -  |  arctan(z/zR)  +  ♦Q,  (la) 

R  =  z(  1  +  (zR/z)2}  (lb) 


is  the  radius  of  curvature  of  the  wavefronts, 


(i)W 

o 

ZR  =  2c 


is  the  Rayleigh  range,  is  a  constant, 


(lc) 


Ac  =  Ac 

s  s 


w 

_0 _ 

w(z) 


1/2 


exp 


(-y2/w2(z)J, 


(Id) 


A  is  a  constant,  and  w(z),  the  spot  size  at  z,  is  related  to  the  waist  w 
s  o 

at  z  =  0  by 


2 


v(z) 


(le) 


It  should  be  remarked  that  the  amplitude  and  phase  of  the  optical  field  are 
assumed  to  be  constant  in  time. 

3 

Following  the  standard  procedure,  the  equations  of  motion  in  terms  of 
the  independent  variable  z  are  deriveable  from  the  Hamiltonian  function  - 
P2(y»Py;t,-E;z),  which  is  given  by 

"z  -  I-  “s-l1  +  <V"C)2  *  ax)’  (2> 

where  E,  the  total  energy  of  an  electron  of  rest-mass  m  and  charge  -|e|,  is 
canonically  conjugate  to  the  time  t,  p^  is  the  momentum  conjugate  to  the 
coordinate  y,  and  a^  is  the  normalized,  total  vector  potential: 

ax  =  awcosh(kwy)sin(kvz)  +  agcos(^  z  -  cot  +  «>),  (3) 

2 

with  a.,  =  |e|A  /(me  ).  Since  the  vector  potential  is  not  an  explicit 

w ,  s  1  1  w ,  s 

function  of  the  coordinate  x,  the  associated  canonical  momentum  is  a 

constant.  Equation  (2)  has  been  derived  by  choosing  Px  =  0,  and  omitting 

2  -2 

terms  on  the  order  of  (me  /E)  and  higher  in  the  expansion  of  the 
Hamiltonian. 

Upon  squaring  the  expression  in  Eq.  (3),  averaging  over  the  wiggler 
period  and  retaining  only  the  slowly-varying  contribution  to  the 
ponderomotive  potential,  it  is  found  that 

p  =  -  -  m yf;  / 1  +  (p  /mc)^  +  i  a^cosh^(k  y)  +  a  a  cosh(k  y) 
c  2E  |  v*y  2  v  '  vJ '  w  s  '  vJ 

x  sin  k-  +  k  )z  -  wt  +  ♦ 

L  c  w' 
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Hamilton's  equations  are  then  given  by 


iL  _ 


dz 


=  v, 


dv  _  .2 

dz  ■ 


sinh(2kwy) 


2k 


af  ■  k«  -  rr[l  + 1 4cosh2ov> ■*  (w)2], 

Zy  c 


3?  =  5^  awasf gCoshCk^ )  cos  ( +  ♦), 


(4a) 

(4b) 

(4c) 

(4d) 


where  y  =  E/(mc  ),  *  =  (<o/c  +  ky)z  -  tot,  v  =  py/(Yinc),  fg  =  JQ(0  -  J^O 
is  the  usual  difference  of  Bessel  functions,  £  =  (ay/2)  /(I  +  a^/2),  and 


k„  = 


a  k 
w  w 

/2  Y 


(5) 


is  the  betatron  wavenumber.  In  arriving  at  Eq.  (4b)  it  has  been  assumed 
that 


ra. 


“i  V 


9  »  2 

1  +  a  /2  a  (k  w  )  , 
w  w '  w  o 7 


«  1, 


(6) 


where  wq  is  the  waist  of  the  optical  field,  Eq.  (le).  Additionally,  the 

2 

wave  term,  proportional  to  a  /y  ,  has  been  dropped  from  Eq.  (4c). 
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III.  Analysis  of  synchrotron-betatron  coupling 

Variation  of  the  optical  phase  $  across  the  wavefronts  underlies  the 

instability  investigated  in  Ref.  1.  Specifically,  the  first  term  in  (la) 

leads  to  a  linear  increase  with  z  of  the  synchrotron  oscillation  amplitude 

of  deeply  trapped  electrons.  Assuming  the  radius  of  curvature  R  of  the 

wavefronts  to  be  constant,  the  growth  rate  of  this  instability  is 
2 

proportional  to  y_/R,  where  yQ  is  the  amplitude  of  the  betatron 
P  P 


2 

oscillation.  In  what  follows  it  is  assumed  that  y g/R  -»  0.  This  is  the 
case  for  electrons  confined  to  the  central  portion  of  the  electron  beam,  or 
in  the  region  where  the  curvature  of  the  wavefronts  is  small.  Assuming 
further  that  the  2-variation  of  ♦  is  small  on  the  scale  length  of  the 
instability  to  be  discussed,  for  a  synchronous  electron  Eqs.  (4c)  and  (4d) 
imply 


Y  = 

Y  =  constant, 

(7a) 

*  + 

♦  =  (n  +  |)n, 

n  =  0 ,  +1,  +2,  ... 

(7b) 

whence,  for  |kwy|  «  1,  Eqs.  (4a)  and  (4b)  integrate  to 

y  =  ygcos(kgz)  +  k^VgSintkgZ) ,  (8) 

where,  from  (5),  kg  =  avkv//2y  is  the  betatron  wavenumber  for  a  synchronous 
electron,  and  y^,  Vg  are  constants.  Upon  substituting  (8)  into  Eq.  (4c) 
one  obtains 


\  aJcosh2(kwy)  +  (yv)2  *  j  a2[l  +  (kpf)2]  +  (wV 

1  2  “2  (.*2  2  2^ 

■*  2  aw  +  Y  IVp  +  veJ’ 


which  is  a  constant  in  z,  whence  the  synchronous  energy  follows  from  Eq. 
(4c)  by  setting  d^/dz  =  0: 

.21 


1  + 


HIl*  <w 


2ck  /« 
w 


(9) 


'e 


For  an  electron  in  the  vicinity  of  the  synchronous  particle, 

Y  =  y  +  6 y»  'J'  =  and  Eqs.  (4c)  and  (4d)  imply  synchrotron 

oscillations  of  Sy  and  of  S\p.  Through  the  Y-dependence  of  the  betatron 
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wavenumber,  Eq.  (5),  it  follows  from  Eqs.  (4a)  and  (4b)  that  under  an 
appropriate  resonance  condition  a  parametric  instability  of  the  coupled 
betatron  and  synchrotron  oscillations  is  then  possible. 

To  examine  this,  the  expression  in  (8)  is  generalized  to 

y  =  yp(z)cos  [(kg  +  e)zj  +  (kg  +  e)-1Vg(z)sin  £(kg  +  e)z],  (10) 

where  ya  and  vQ  are  now  functions  of  z  and  e  is  a  small  wavenumber  shift  to 
P  P 

be  determined.  Assuming  |dyg/dz  |  ~  ey^,  |dVg/dz |  ~  £Vg,  and  omitting  terms 
2 

on  the  order  of  e  ,  substitution  of  (10)  into  Eqs.  (4a)  and  (4b)  yields 
(dVg/dz  -  ek^y^cos  |(kg  +  £)zj  -  (kgdyg/dz  +  evg)sin[(kg  +  e)zj 

=  k^  |ygCos  [(kg  +  e)zj  +  (kg  +  e)-1VgSin  [(kg  +  e)zj|(8y/r)- 


To  obtain  an  expression  for  &y(2)  to  be  used  in  Eq.  (11),  Eqs.  (4c) 
and  (4d)  are  perturbed  about  the  synchronous  values  (7a)  and  (7b): 

fes^ky  «  [l  +  ia2(1  +  ky)  +  ;2v2] 


-%[i  +  jaJ(l  +  kjy2)]«  y, 


(12a) 


,  ua  a  f. 

d  «.  ,  ,4n  w  s  B 

31  Y  ‘  l'1>  ~ 


(i  *  \  kjy2)s*, 


(12b) 


where  n  denotes  an  integer  from  the  set  defined  in  (7b).  Substituting  (10) 
and  the  forms 


Sy  =  Ay(z)cos  j^2(kg  4  e)z  4  h] » 
=  AY(z)sin  [2(kg  +  e)z  +  hj  > 


(13a) 


(13b) 
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where  n  is  a  constant,  into  (12)  and  assuming  that 


i  *-1/2 

Vs «  T  ’ 


(14) 


one  obtains  an  identity  whose  consistency  requires  that  Ay  and  be 


independent  of  z 

and  tha  t 

2(ke  + 

e)  =  k  , 
syn’ 

(15) 

where 

k 

syn 

2 

m  a  a  fD  , 
~  w  s  B  (. 

L  ,'4  2  { 

2y  c 

i  1  2)11/2 

1  *  2  avJJ 

is  the  synchrotron  wavenumber.  In  addition,  it  is  required  that  the 
integer  n  in  Eq.  (12b)  be  odd.  It  should  be  noted  that  the  constraint 
indicated  by  (14)  is  consistent  with  the  neglect  of  the  radiation  field 
term  in  Eq.  (4c)  and  the  parameter  restriction  in  (6). 

Now  substituting  (13a)  into  Eq.  (11)  one  again  obtains  an  identity. 
Assuming  y^,  v^  -  exp(Az),  the  consistency  condition  reduces  to 


A  = 


(16) 


which  is  the  growth  rate  of  the  instability. 


IV.  Discussion 

It  is  clear  that  the  source  of  the  instability  discussed  here  lies  in 
the  energy  dependence  of  the  betatron  wavenumber,  Eq.  (5).  Previous 
studies  of  the  synchrotron-betatron  coupling  have  neglected  this 
dependence.  The  results  of  Refs.  1  and  2  indicate  a  linear  growth  in  the 
amplitude  of  the  synchrotron  motion  (Ay,  Avj/)  for  deeply  trapped  electrons. 
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On  the  other  hand,  for  the  parametric  instability  discussed  herein,  Ay  and 


A<|/  are  constants,  whereas  and  v^  grow  exponentially  fast.  In  all  three 

cases,  however,  the  resonance  condition  is  essentially  k  =  2kfl. 

syn  p 

From  Eqs.  (15)  and  (16)  it  follows  that  for  a  sufficiently  large 
wavenumber  mismatch  e  the  growth  rate  vanishes.  Maximum  growth  is  obtained 
for  e  =  0: 

Xmax  ■=  k  k0|(Ay/y)cosT1|.  (17) 

Since,  at  the  order  of  approximation  employed  here,  Ay  and  Ai|/  are 
constants,  the  parametric  instability  is  expected  to  affect  the  extraction 
efficiency  of  the  FEL  mainly  through  ir  '.reased  electron  beam  radius  and 
particle  detrapping  from  the  ponderomotive  potential,  such  as  that 
described  in  Refs.  1  and  2,  would  result  once  the  betatron  amplitude  is 
sufficiently  large.  It  must  be  remarked,  however,  that  if  the  variation  of 
the  optical  field  amplitude  or  phase  over  the  scale  length  indicated  by 
(17)  is  large  enough,  then  the  resonance  (15)  may  be  crossed  without  an 
appreciable  effect  on  the  efficiency. 
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